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We present our preliminary result for the charmed quark mass, which follows from taking the D s and K meson 
masses from experiment and ro = 0.5 fm (or, equivalently Fk = 160 MeV) to set the scale. For the renormalization 
group invariant quark mass we obtain M c = 1684(64) MeV, which translates to m c (m c ) = 1314(40) (20) (7) MeV 
for the running mass in the MS scheme. Renormalization is treated non-perturbatively, and the continuum limit 
has been taken, so that the only uncontrolled systematic error consists in the use of the quenched approximation. 



1. INTRODUCTION 

Quark masses are among the fundamental pa- 
rameters of the Standard Model, and yet not di- 
rectly accessible to experimental measurement. 
Their determination in terms of experimental 
quantities requires a good quantitative control 
of the non-perturbative dynamics of QCD, which 
may be achieved by numerical simulations of the 
lattice regularized theory. Extending previous 
work on the A-parameter and the strange quark 
mass m quenched QCD g@] we report here on a 
precise determination of the charm quark mass, 
taking the experimental value of the D s meson 
mass as essential new input. For previous work 
on the charm quark mass we refer to two recent 
papers || and the references given there. 

2. STRATEGY 

We use O(a) improved Wilson quarks with the 
Wilson gauge action. All necessary renormaliza- 
tion constants and improvement coefficients are 
known with non-perturbative accuracy in the rel- 
evant range of bare couplings. We refer to |Q] for 
a review and further references. 

2.1. Setting the parameters 

In quenched QCD it is convenient to measure 
dimensionful quantities in units of the scale ro Q , 
since ro/a has been determined very precisely for 
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a large range of bare couplings (^J^]. To ob- 
tain the scale in physical units we always use 
ro = 0.5 fm, which is roughly equivalent to set- 
ting the scale with F K = 160 MeV @. The renor- 
malization group invariant (RGI) strange quark 
mass has been obtained in , 

r M s = 0.348(13) =4> M s = 138(6) MeV, (1) 

taking mx from experiment and the mass ratio 

M s /M = 24.5 ± 1.5, M = \{M U + M d ), (2) 

from chiral perturbation theory ]8j . Up to intrin- 
sic 0(a 2 ) ambiguities, the connection between M s 
and the bare strange quark mass can be estab- 
lished by combining the results of refs. [|l]|j], so 
that no tuning is required to satisfy the renor- 
malization condition (^). The bare charm quark 
mass may then be fixed by matching directly the 
experimental result mrj s = 1969 MeV (l(J, which 
translates to romo, — 4.98. This is justified since 
electroweak effects on D meson masses are ex- 
pected to be small on the scale of the expected 
statistical errors. 

2.2. Some technical details 

The simulations were done at four values of 0, 
6.0 < f3 < 6.45, which correspond to lattice spac- 
ings in the range 0.05 — 0.1 fm. Spatial lattice vol- 
umes ranged from 16 3 to 32 3 , and the Euclidean 
time extent from 40 to 80. The linear extent of 
the spatial volume was never less than 1.5 fm in 
physical units, so that finite volume effects on the 
D s meson mass can be safely neglected. Quark 
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propagators were computed for two mass parame- 
ters in the strange quark region, and three values 
around the expected charm quark mass. To com- 
pute meson masses we used correlation functions 



derived from the Schrodinger functional (SF) |11 



In the pseudoscalar channel these are given by 



fx(x ) 



i« 6 



^(c(x)r xS (x)C s (y) 7 5Cc( Z )),(3) 



where r x stands for 7075, 75 for X = A,P re- 
spectively. The pseudoscalar meson mass is then 
obtained by looking for a plateau in the effective 
masses, 

amf s (x + |a) = In {f x (x ) / fc(x + a)} , (4) 

as a function of Euclidean time Xq. Autocorre- 
lation times were found to be completely negligi- 
ble, so that the produced 0(100) configurations 
at each (3- value are statistically independent. An 
example for an effective mass plot is given in fig. 1. 
Excited states are expected to contribute signifi- 
cantly when either xq or T — xo becomes small. 
It is possible to estimate these contributions and 
thereby quantify the resulting deviation from the 
true plateau value. Requiring this effect to be 
smaller than 0.5%, the plateau region in fig. 1 
was identified with the interval 25 < xo/a < 35; 
the meson mass was then estimated by averaging 
the effective mass in this interval. 



2.3. Bare and renormalized quark masses 

Given the pseudoscalar masses at the chosen 
parameters, an interpolation in some bare charm 
quark mass is required to match the experimen- 
tal D s meson mass. Besides the bare mass too 
we also consider two alternative definitions of the 
bare charm quark mass, which derive from the 
PCAC relation. First we use the heavy-light ax- 
ial current and density to define, 



dofAjxp) + c A adQdofp(xo) 
2fp(x ) 



(5) 



which upon renormalization yields the average 
charm and strange quark mass. With our param- 
eter choices yet another bare charm quark mass 
can be obtained from the PCAC relation involv- 
ing a hypothetical strange quark, assumed to be 
mass degenerate with the charm quark. We de- 
note this PCAC mass as m cc but one should bear 
in mind that its definition does not involve flavour 
singlet operators. 

The renormalized and 0(a) improved charm 
quark mass tor can now be obtained in various 
ways. Choosing the SF scheme we first obtain 



tor = Zp 1 Zm q (l + b m am q ), 



(6) 



where m q is the subtracted bare charm quark 
mass. From the PCAC mass to cs , one gets 

«ir+tor :S = Zp X Z K {1 + (6 A - 6p)aTO q } 2to cs , (7) 

where m q denotes the average of the bare sub- 
tracted charm and strange quark masses, and we 
have assumed xq — T/2, which is kept approx- 
imately constant in physical units. Finally, we 
consider 



Z p 1 Z A {1 + (6 A - b P )am q } m CCl 



(8) 



and convert to the RGI quark mass using the 
flavour-independent ratio 0] 



M/tor = 1.157(15). 



(9) 
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Figure 1. Effective pseudoscalar mass in lattice 
units vs. Xo/a at ft = 6.1. 



3. RESULTS AND CONCLUSIONS 

After subtraction of the strange quark mass (|l|) 
in eq. ^ we have three definitions of r$M c which 
should coincide up to terms of 0(a 2 ). In each case 
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we extrapolate to the continuum with an ansatz 
of the form 



ALPHA collaboration 



r M c = c + ci(a/r ) 2 



(10) 



excluding the data point at the coarsest lattice 
spacing ((3 = 6.0). As can be seen in fig. 2, the 
results are nicely compatible with each other in 
the continuum limit, although the differences at 
finite cutoff are quite large. As our best (prelim- 
inary) estimate we quote 

r M c = 4.26(16) =*> M c = 1 684(64) MeV. (If) 

Using the four-loop anomalous dimension and (3- 
function in the MS scheme |L2|], and r A^ 
0.586(48) §, we obtain 



""MS 



m c (m c ) = 1314(40) (20) (7) MeV, 



(12) 



where the first error stems from the RGI quark 
mass, the second is induced by the error on A-^g, 
while the last number is the difference between 
three- and four-loop perturbative evolution down 
to the charm scale. 

We conclude that precise results can be ob- 
tained for charmed observables using the same 
techniques as for the light quarks, provided that 
the continuum extrapolation is performed. We 
emphasize that we have used the quenched ap- 
proximation, which is known to lead to inconsis- 
tencies in the hadron spectrum at the f 0% level. 
However, the present result constitutes a signifi- 
cant improvement over previous lattice estimates 
and will be a useful reference point for more real- 
istic calculations involving dynamical quarks. For 
the time being we speculate that the mass ratio 
M c /M s = 4.26(f6)/0.348(13) w 12.2±1.0may be 
less sensitive to quenched scale ambiguities than 
the masses themselves. 
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Figure 2. Continuum extrapolation of the RGI 
charm quark masses. 
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